Magnetic Moment of the Fragmentation- Aligned 61 Fe(9/2 + ) Isomer. 
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We report on the g factor measurement of the isomer in 61 Fe (E* = 861 keV). The isomer was 
produced and spin-aligned via a projectile-fragmentation reaction at intermediate energy, the Time 
Dependent Perturbed Angular Distribution (TDPAD) method being used for the measurement of 
the g factor. For the first time, due to significant improvements of the experimental technique, an 
appreciable residual alignment of the isomer has been observed, allowing a precise determination 
of its g factor: g — —0.229(2). Comparison of the experimental g factor with shell- model and 
mean field calculations confirms the 9/2 + spin and parity assignments and suggests the onset of 
deformation due to the intrusion of Nilsson orbitals emerging from the vgg/2- 

PACS numbers: 21.10.Ky, 21.60.-n, 23.20.En, 25.70.Mn 



The measurement of electromagnetic (EM) moments 
has traditionally played a central role in the critical evalu- 
ation of nuclear structure models since they elucidate the 
single-particle nature (magnetic moments) and the shape 
(quadrupole moments) of the nuclear state under inves- 
tigation. Their particular behavior around nuclear shell 
closures, approaching Schmidt values for magnetic mo- 
ments and small values for quadrupole moments, makes 
them good tools to investigate shell closure near and far 
from P stability. 

One of the main restrictions in the study of EM mo- 
ments is the necessity to obtain oriented isomeric ensem- 
bles. In particular, for the study of neutron rich nu- 
clei, that are mostly produced via fragmentation reac- 
tions at high or intermediary energies, the mechanism to 
produce oriented states is not yet well understood. The 
first observation of spin alignment of isomeric states in 
a projectile-fragmentation reaction at an energy of 500 
MeV/u was reported by Schmidt-Ott et al 0. At pro- 
jectile energies below 100 MeV/u, so called intermediate 
energies, spin alignment of isomeric beams was reported 
recently by Georgiev et al 2] in the study of isomeric 
states in the 68 Ni region. However, only a very small 
residual alignment observed in the decay of the 67 Ni and 
the 69 Cu isomers was reported. 

In this letter, we report an important experimental 
achievement in the study of EM moments of neutron rich 
isomers. For the first time, high quality data could be 
obtained for an isomer in the neutron rich 61 Fe nucleus, 
located near the controversial N=40 subshell closure la, 



We focus here on the role played by the vgg/2 orbital 
in the low-energy level structure of nuclei near N = 40 in 
the particular case where this orbital manifests itself as 
an isomeric state. The 61m Fe, with 35 neutrons, is one 
of the lightest nuclei exhibiting such an isomeric state 
Pi/2 = 250(10) ns) at a rather low excitation energy 
(E* = 861 keV). The 9/2+ tentatively assigned spin and 
parity is based on systematics @. The measurement of 
the g factor of the isomeric state in 61 Fe provides in- 
formation about its structure and can also confirm the 
suggested spin/parity. 

The nuclei of interest were produced following the frag- 
mentation of a 54.7 MeV/u 64 Ni beam accelerated at the 
GANIL facility, with a mean intensity of 7 • 10 11 pps, 
impinging on a 97.6 mg/cm 2 thick 9 Be target placed 
at the entrance of the LISE spectrometer [lfj. In or- 
der to decrease the in-flight decay of the isomer, the de- 
tection setup (Fig. was positioned at the first focal 
plane of the LISE spectrometer (time of flight « 200 
ns). A 300 [im thick removable silicon detector was used 
to optimize the selection of 61 Fe fragments by means of 
energy-loss vs time-of-flight identification. Once the se- 
lection was performed, a 50 /im thick plastic scintillator, 
placed in front of the catcher foil, was used to provide 
the t = signal for the subsequent TDPAD measure- 
ment. The choice of a thin plastic scintillator instead of 
a silicon detector was made for two reasons. First, the 
fragments are selected fully stripped with the LISE spec- 
trometer and, in order to preserve the initial alignment 
until the eventual stop in a catcher foil, it is important 
to avoid as much as possible the capture of electrons that 
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FIG. 1: Schematic drawing of the TDPAD experimental set- 
up. The beam passes through a 50 /im plastic scintillator 
before being stopped in a 500 /im Cu foil. 



can completely destroy the orientation of the nuclear en- 
semble. The electron pick-up upon pas sing through the 
scintillator is estimated to be below 2% |ll|, much lower 
than that in a 300 fim thick silicon detector (60-70%) [2j. 
Secondly, the accepted secondary beam rate (between 17 
and 80 kHz for the present experiment) can be one order 
of magnitude higher than with a silicon detector. As a 
catcher for the reaction products we used an annealed 
high-purity 500 \xm thick Cu foil. Iron ions have the 
same electronegativity and similar atomic radius as Cu 
atoms and hence Cu, with its cubic lattice structure, is 
expected to provide a perturbation-free environment for 
implanted Fe fragments. The Cu foil was placed between 
the poles of an electromagnet that provided a constant 
magnetic field ~$ in the vertical direction. 

The Larmor precession of the initially aligned isomeric 
spins in the applied field of about 0.7 T, was monitored 
with four coaxial Ge detectors placed in the horizontal 
plane around the Cu foil as shown in Fig. Time spec- 
tra were collected, having as start the signal due to the 
ion passage through the plastic scintillator and as stop 
the signal given by the detection of a prompt or delayed 
7-ray. To extract the precession pattern out of the time 
spectra, data from detectors positioned at 90° with re- 
spect to each other was combined to generate the R(t)- 
function: 



R(t) 



I 12 (8,t)-el u (6+i,t) 
I 12 (9,t)+el 34 (e+i,t) 
A 2 B%{t = 0)cos{2(u L t + a- 9)). 



(1) 



with o?i=-ffi^ and £*=-§ (1 - f§) both depending 
on the g factor. Ii 2 and I34 are the summed intensities of 
detectors placed at 180 degrees (see figure P); A 2 is the 
radiation parameter of the 7-ray transition; B 2 is the sec- 
ond component of the orientation tensor describing the 
initial orientation of the selected isomeric ensemble; 9 
equals 7r/4; e is the relative efficiency between the four 
respective detectors in the present setup, and a is the ro- 
tation angle of the aligned ensemble symmetry axis with 
respect to the beam axis induced when passing through 



the two dipole magnets of LISE spectrometer 0. The 
data acquisition was validated on an event by event ba- 
sis by the coincidence between a heavy ion signal from 
the plastic scintillator and a delayed 7, registered by one 
of the germanium detectors, within a time window of 3 
fis. In order to diminish the number of accidental co- 
incidences in this time window, we have used a pack- 
age suppresser that brought on target 1 out of 10 ion 
packages provided by the accelerator with a frequency of 
10.5 MHz. As a consequence, the detected alignment in- 
creased by a factor of ~ 3 compared to a measurement 
without package suppression. 

Certain physical effects need to be taken into account 
when deducing a precise value for the g factor from the 
observed Larmor frequency. The distribution of the mag- 
netic field over the beam spot, the paramagnetic amplifi- 
cation of the applied magnetic field and the Knight shift, 
all can induce minor modifications of the Larmor fre- 
quency. In order to avoid systematic errors due to these 
corrections and to validate the experimental setup and 
method, we have measured the Larmor precession under 
identical conditions for isomers in two Iron isotopes: the 
j7r _ 20+, E* — 6527 keV isomer in 54 Fe having a known 
g factor, g(10+) = +0.7281(10) H2, and the 7 7r =(9/2+), 
E* = 861 keV isomer in 61 Fe decaying via a (M2) transi- 
tion of 654 keV in cascade with a (Ml) transition of 207 
keV to the ground state (see fig. EJ). The effective value 
of the magnetic field extracted from the R(t) function 
of 54m Fe (eq. HJ) is 6800(40) Gs and it includes all the 
corrections mentioned above. 

54m Fe, decaying by a stretched E2-cascade, served 
also as a probe to measure the produced alignment of 
the isomeric ensembles as a function of their angular 
momentum distribution |l3j| . For the selection of the 
fragments in the wing of the momentum distribution 
ip fragment < Pprojectiie) , the amplitude of the R(t) func- 
tion (fig. EJ) yields a large negative alignment, —12.5(9)%, 
the sign being in agreement with the predictions of a 
kinematical fragmentation model [lj, [l5j . 

For 61m Fe, we have studied the alignment for two en- 
sembles selected, respectively, in the center and in the 
outer wing of the longitudinal momentum distribution. 
The alignment of the isomeric state, deduced from the 
amplitude of the R(t) functions of the 654 keV decay is 
+6.2(7)% for the central selection and -15.9(8)% for the 
wing selection, assuming a pure M2 transition. 

The R(t) functions of the 207 keV and the 654 keV 
7-rays de-exciting the 61m Fe isomer exhibit an opposite 
phase and the ratio between their amplitudes is 1.43(16). 
Using realistic GEANT simulations 0] and assuming 
that the 654 keV and 207 keV transitions have pure M2 
and Ml multipolarities, respectively, and the level se- 
quence is 9/2+^ 5/2~^ 3/2 _ , we estimated this ratio to 
be 1.30(6). A spatial distribution close to the fragmented 
beam spot of the emitting 61 Fe isomers is assumed as 
well in order to take into account the geometrical factor 
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FIG. 2: The R(t) functions for: a.- 61 Fe(center, 654 keV); b.- 
61 fe(wing, 654 keV); c.- 61 J Fe(wing, 207 keV); d.- 54 i^e. The 
R(t) function for 207 keV is consistent with that for the 654 
keV but with opposite phase, as expected from the assumed 
multipolarities of the corresponding 7 transitions. 



of the detection setup. The good agreement between the 
two ratios (experimental and simulated) indicates that 
the multipolarities and sequence of the 7 transitions used 
for the determination of the experimental alignment are 
correct. 

The measured half-life of the 61 Fe isomer, Ti/ 2 = 
245(5) ns, is in good agreement with the previous mea- 
surements 9J. For the g factor, a value of —0.229(2) 
was extracted by a % 2 fit of the R(t) functions using ex- 
pression Q with both frequency and initial phase, a, 
depending explicitly on g. The error on the fitted g fac- 
tor includes the errors on the effective field. The fitted 
R(t) functions were constructed using the 654 keV di- 
rect decay of the isomeric state for the central and wing 
selection of the fragments momentum, respectively (fig. 

In Fig. |3| we present the measured g factors of 9/2+ 
isomeric states around N = 40, including also the g fac- 
tor of 61m Fe. The comparison with the other g factors 
of known 9/2+ isomeric states in the region strongly sup- 
ports the 9/2+ spin and parity assignment for the 61m Fe 
isomer. One can observe the symmetry with respect to 
Z = 28 of the g factors for N = 35 chain and the increase 
of g factor values when one goes away from Z = 28, in- 
dicating an increase of core polarization effects. 

In order to extract detailed information about the va- 
lence orbital occupation for the isomeric state, we have 
performed large scale shell model (LSSM) and mean-field 
calculations using the Hartree-Fock-Bogoliubov (HFB) 
formalism with an effective D1S Gogny nucleon-nucleon 
force 0. 

The shell model calculations were performed using the 
ANTOINE code [13] of the Strasbourg group in a valence 
space composed of f,p and gg / 2 active orbitals and having 
a closed 48 Ca core. The interaction used is described in 
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FIG. 3: g factor systematics around N=40 for known 9/2 + 
isomeric states. The data are taken from 0, lr?l Hill . The 
unusual value for the g factor of 67 Ni could be explained by 
proton excitations across Z—2S. 



[8j and in the references therein. 

The calculated free g factor is 0(9/2+) = -0.277. The 
effective value is g(9/2 + ) = —0.1828 if a quenching fac- 
tor of 0.7 is used for the nucleon spin g factor. The 
value of the quenching factor is quite arbitrary because 
at present there is no systematical comparison between 
experimental data and shell model calculations into the 
considered space. Recently, it was shown that for the sd 
and pf shells, the configuration mixing within the shells 
is enough to fully account for the observed magnetic mo- 
ment, without the use of t he q uenching factor needed for 
Gamow- Teller beta decay |21|. The 9/2+ state is calcu- 
lated to be at 720 keV and its wave function is a mixture 
of a large number of configurations, but the mean occu- 
pation of the vgg/2 orbital is ~ 1. The good agreement 
between the calculated and the experimental value of the 
g factor indicates that the isomer is a 9/2+ state, gen- 
erated by a neutron in the vgg/2 orbital whilst its very 
mixed wave function could be related to a deformed po- 
tential. 

For a better understanding of the variation of the sin- 
gle particle orbitals near the Fermi surface as a func- 
tion of the nuclear deformation, we have calculated the 
ground-state and excited states within the Hartree-Fock- 
Bogoliubov (HFB) formalism. The calculation beyond 
the mean-field approximation was performed using the 
Generator Coordinate Method (GCM) with the Gaussian 
Overlap Approximation [2^ . The ground state is found 
to be a K 1 * = 3/2~ prolate deformed state, fa — 0.187, 
whereas an excited state K* — 9/2+ is predicted with 
an excitation energy E* = 1640 keV and an oblate defor- 
mation, fa = —0.227. Fig. 0] presents the HFB energies 
of the states located near the Fermi level in 61 Fe as a 
function of the deformation parameter, fa. 

The free g factor of the J v = 9/2+ state built from the 
K* = 9/2+ HFB state is found to be g(9/2+) = -0.3477 
(without a quenching factor) and it remains practically 
unchanged after applying the GCM method. 

The comparison between the experimental value and 
the theoretical predictions, LSSM and HFB, indicates 
a deformed potential for the 9/2+ isomeric state. One 
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FIG. 4: Potential energy curves of the ground state and of 
the first three excited states of 61 Fe as functions of the de- 
formation parameter, 0%. 



order to understand the coupling of the unpaired parti- 
cle to the 60 Fe core, it is important to measure the sign of 
the quadrupole moment of the isomeric state, requiring 
a spin polarized isomeric beam |25l l2fij | . 

The appreciably large residual alignment measured for 
the 61 Fe and 54 Fe fragments indicates that fragmen- 
tation reactions at intermediate energies can provide a 
powerful tool to align ensembles of nuclear isomers of yet 
more exotic nuclei, thus facilitating the determination of 
EM moments in neutron rich nuclei and allowing the in- 
vestigation of nuclear structure away from stability. 
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can notice that the two extreme HF orbitals (9/2+ and 
1/2+) generating from the spherical vg§/2 orbital (fig. |3J 
are close to the Fermi surface at moderate deformations 
but with opposite sign. In the one particle plus rotor 
model |23| . this particularity provides two possible ways 
of creating a J n = 9/2 + state, depending on how the 
unpaired neutron is coupled to the rotating 60 Fe core. 
In the strong coupling scheme [24| , the unpaired neutron 
spin couples to the deformation of the core and the 9/2 + 
state has the projection on the symmetry axis K = 9/2. 
The second possibility is that the unpaired neutron spin 
couples to the rotation of the core and the total spin 
projection on the symmetry axis is degenerated between 
K = 1/2 and K = —1/2 due to the Coriolis interaction, 
so K is no longer a good quantum number. Instead, the 
projection on the rotation axis becomes a good quan- 
tum number and the spin of the physical state can be 
9/2+ . We have then calculated the low energy states in 
the frame of this model for 61 Fe. Taking the deforma- 
tions indicated by the HFB calculations, for the two cases 
presented above, the lowest state having a positive parity 
has a 9/2 spin and an excitation energy of about 850 keV. 
The g factors are similar, which was predictable because 
the g factor should not be sensitive to the deformation. 

In conclusion, the measured g factor is in very good 
agreement with the assigned 9/2 + spin and parity. From 
the comparison with LSSM and HFB calculations there 
are indications that this state is characterized by a de- 
formed potential, induced by a steep lowering of the Nils- 
son orbitals emerging from the spherical vgg/2 neutron 
orbital as a function of deformation parameter, 02- In 
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